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ABSTRACT 

We present new spectroscopic observations of the gravitational arcs and the brightest cluster galaxy 
(BCG) in the cluster MS2137-23 (z = 0.313) obtained with the Echellc Spectrograph and Imager on 
the Keck II telescope. We find that the tangential and radial arcs arise from sources at almost identical 
redshifts [z — 1.501,1.502). We combine the measured stellar velocity dispersion profile of the BCG 
with a lensing analysis to constrain the distribution of dark and stellar matter in the central 100 kpc of 
the cluster. Our data indicate a remarkably flat inner slope for the dark matter profile, pd oc r _/3 , with 
/3 < 0.9 at 99% CL. Steep inner slopes obtained in cold dark matter cosmological simulations - such as 
Navarro Frenk & White ((3 = 1) or Moore (1.5) universal dark matter profiles - are ruled out at better 
than 99%CL. As baryon collapse is likely to have steepened the dark matter profile from its original 
form, our data provides a powerful test of the cold dark matter paradigm at the cluster mass scale. 



Subject headings: gravitational lensing — 
— galaxies: formation 



galaxies: elliptical and lenticular, cD 
— galaxies: structure 



galaxies: evolution 



1. INTRODUCTION 



A fundamental result arising from cold dark matter 
(CDM) numerical simulations is that the density profiles 
of DM halos are universal in form across a wide range 
of mass scales from dwarf galaxies to clusters of galaxies 
(Navarro, Frenk & White 1997, hereafter NFW). Internal 
to some scale radius r sc , the dark matter profile assumes 



a power law form, pd oc 



Whilst there is some dis- 



pute amongst the simulators about the precise value of (3 
with values ranging from 1.0 to 1.5, (Moore et al. 1998, 
hereafter M98; Ghigna et al. 2000, Power et al. 2002), a 
clear measurement of (3 in a range of objects would offer a 
powerful test of the CDM paradigm. 

The largest observational effort in this respect to date 
has been via dynamical studies of low surface brightness 
(LSB) and dwarf galaxies, suggesting softer (f3 < 1) DM 
cores than expected on the basis of the numerical simu- 
lations (e.g. de Blok & Bosma 2002, Salucci & Burkert 

2000) , although the issue remains somewhat controversial 
(e.g. van den Bosch & Swaters 2001). Similar tests have 
recently been extended to regular spiral (Jimenez, Verde & 
Oh 2002) and elliptical galaxies (Treu & Koopmans 2002). 
Some observational constraints are available at the scale 
of massive clusters, from lensing (e.g. Tyson, Kochanski 
& Dell' Antonio 1998; Williams et al. 1999; Smith et al. 

2001) , X-ray analysis (Mahdavi & Geller 2001) and dy- 
namics of cD galaxies (Kelson et al. 2002). Since massive 
clusters probe a totally different scale and physical condi- 
tions than galaxies, it is crucial to understand their mass 
distribution to test the universality of the DM profiles. 

In this paper we present the first application of a new 
method to determine the luminous and dark mass distri- 
bution in the inner regions of massive clusters with giant 



arcs around a central BCG. The method combines lens- 
ing analysis with stellar kincmatical measurements of the 
BCG. The two ingredients provide complementary infor- 
mation on the relevant scales (~ 100 kpc), allowing us to 
disentangle the luminous and dark components of the total 
mass distribution. 

We have chosen the cluster MS2137-23 as a first appli- 
cation of our method since it is an approximately round 
system, has an isolated BCG and a very well-studied arc 
system. Fort et al. (1992) first pointed out the potential 
significance of the radial and tangential gravitational arcs 
as a means of constraining the mass distribution on ~100 
kpc scales, and mass models have been developed subse- 
quently (Mellier et al. 1993, Hammer et al. 1997; hereafter 
M93, H97). The redshifts of the radial and tangential arcs 
were predicted to lie in the range 1< z <2 (M93). A key 
issue in the earlier work is whether the radial arc is, in fact, 
a lensed feature. M93 and Miralde-Escude (1995; hereafter 
ME95) also point out the importance of determining the 
stellar velocity dispersion profile of the BCG to weigh the 
stellar contribution to the mass. 

Following the earlier suggestions, we present new ob- 
servations of the cluster MS2137-23 with the Keck II 
telescope. We provide spectroscopic confirmation of the 
arcs and measure a velocity dispersion profile for the cen- 
tral BCG. The spectroscopic data are used together with 
archival HST images to constrain the luminous and DM 
distribution of the cluster. In the following, r is the radial 
coordinate in 3-D space, while R is the radial coordinate 
in 2-D projected space. We adopt Ho=65 km s _1 ,Mpc -1 , 
fi TO = 0.3 and = 0.7 for the cosmological parameters. 



1 Based on observations collected at the Keck Observatory, which is operated jointly by the California Institute of Technology and the University 
of California, and with the NASA/ESA HST, obtained at STScI, which is operated by AURA, under NASA contract NAS5-26555. 
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2. OBSERVATIONS 

2.1. i-fecfc Spectroscopy 

We observed MS2137-23 using the Echelle Spectrograph 
and Imager (ESI; Sheinis ct al. 2002) on the W. M. Keck- 
II telescope for a total integration time of 4900s (2 x 1800s 
+ 1300s) on 21 July, 2001. The seeing was 0"6 and the 
1725 x 20" slit was oriented North-South to include the 
BCG, radial arc, and tangential arc (Figure 1). The spec- 
troscopic goals were two-fold: a determination of the red- 
shift of the arcs and a measurement of the internal kine- 
matics of the central galaxy. An iraf package was devel- 
oped for the specific task of removing echelle distortions 
while preserving the 2-D shape of the spectrum essential 
for the latter goal (EASI2D, Sand et al. 2002, in prepara- 
tion). The instrumental resolution of ESI was measured 
from unblended sky lines to be 30±7 km s _1 . 

The velocity dispersion profile of the BCG (Figure 2) 
was measured using spectral templates based on several 
G-K giants observed with a 0'.'3 slit. These were smoothed 
to match the instrumental resolution of the l'/25 slit and 
redshiftcd to that of the BCG (z=0.313). Analysis was 
restricted to a region around the G band by virtue of the 
high signal/noise and minimal effect of sky line residu- 
als, using the Gauss-Hermite pixel-fitting software (van 
der Marel 1994). The error bars shown in Fig. 2 represent 
a combination of uncertainties arising from Poisson noise 
and systematics, the latter determined from the scatter 
observed using different templates and continuum fits. 




Fig. 1 — HST F702W image of MS2137-23. The rectangular box 
shows the dimensions (1"25 X 20") and position of the ESI slit used 
to make the observations. The BCG, the radial and tangential arc 
are clearly visible at the bottom, center, and upper end of the slit. 

The high spectral resolution of ESI proved crucial in 
clinching the redshifts of the arcs as the emission lines are 
located in a crowded region of OH sky background. The 
two top panels in Fig. 2 show the relevant portion of the 
ESI spectra for the tangential and radial arcs, with the 
observed emission lines identified as the [Oil] doublet at 
z = 1.501 and z = 1.502 respectively. The [Oil] doublet is 
clearly resolved for the tangential arc and it is reasonable 
to suppose that the missing component for the radial arc is 



obscured by sky emission. No other lines are detected on 
either spectra, down to the blue cutoff of ESI at ~ 3900A. 
This makes it unlikely that the single line observed for 
the radial arc is any of the common lines such as Ha, H/3, 
[OIII]4959, 5007, CIV1549, HeII1640, C[III]1909 because 
bluer lines would be detected assuming typical flux ratios. 
The identification of the line with Lya at z — 6.66 is also 
unlikely given that the arc is detected in the HST F702W 
image (see below). Detailed modeling of MS2137-23 based 
on the image configurations predicted that the sources for 
the arcs would be at nearly the same redshift (M93, H97). 
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Fig. 2. — Spectroscopic results: (Top) Strong emission lines de- 
tected in the spectra of the tangential and radial arcs. These are 
identified as [OII]3726,3729 at z = 1.501 and z = 1.502 respectively 
(marked). It is argued that the missing 3729 A line in the radial arc 
is obscured by sky emission. (Bottom) Stellar velocity dispersion 
profile of the brightest cluster galaxy (points with error bars). The 
superimposed histogram shows the profile of the best fitting Jaffe 
+ generalized NFW mass model (see Sec. 3 for details), taking into 
account the effects of seeing (0'.'6), slit width, and radial binning. 



2.2. Hubble Space Telescope Imaging 

HST WFPC2 images of MS2137-23 (GO 5402, PI: 
Gioia), comprising 10 F702W exposures with a total in- 
tegration time of 22.2ks, were used to measure the surface 
photometry of the BCG and to locate arc positions. The 
exposures were reduced in a standard way, using the iraf 
package drizzle (Fruchter & Hook 2002). 

The circularized surface brightness profile (in agreement 
with H97) was obtained using the iraf task ellipse and 
a fit performed as described in Treu et al. (1998, 2001a) 
taking into account the HST point-spread function. The 
best fitting i? 1 / 4 parameters are summarized in Table 1. 
To convert from F702W magnitudes to V magnitudes a 
k-color correction was calculated using the same method 
as Treu et al. (1999). Rest frame photometric quanti- 
ties were corrected for Galactic extinction using Apwiw 
= 2.435E(B-V)=0.122 (Schlegel et al. 1998). 



The dark matter profile of MS2137-23 
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Redshift (BCG) 
Radial critical line 
Tangential critical line 
(l-b/a) e 
F702W (mag) 
SBe,F702W (mag/arcscc 2 ) 

Re,.F702W 

My (mag) 

SB e y (mag/arcsec 2 ) 
R e .v (kpc) 



0.313 ± 0.001 
4"5 ± Of! 3 
15'/35 ± 0'.'20 
0.17 ±0.01 
16.48 ±0.07 
23.58 ±0.34 
5'/02 ± Of! 50 
-24.38 ±0.09 
22.76 ±0.34 
24.80 ± 1.68 



TABLE 1: Relevant spcctro-photomctric quantities 

3. LUMINOUS AND DARK MATTER DISTRIBUTION 

We now combine the observed spectroscopic and pho- 
tometric data to constrain the matter distribution in the 
central region of MS2 137-23. First we introduce a simple 
two-component spherical mass model comprising the stel- 
lar mass of the BCG and a DM halo (Sec. 3.1). We then 
constrain the free parameters of the model using the po- 
sition of the critical lines (3.2) and the velocity dispersion 
profile (3.3). 

3.1. Mass model 
For the luminous component we used a Jaffe (1983) 

M L rj 



(1) 



Tangential and radial critical curves occur when A t = 
and A r = 0, respectively. In practice, the position of the 
tangential arc constrains the total enclosed mass, while 
the position of the radial arc constrains its derivative. A 
proper account of ellipticity is essential for detailed lens 
modeling where the shape, magnification, and morphol- 
ogy of multiple lensed images is being reproduced. ME95, 
however, considered several different simple mass models 
where only the position of the radial and tangential criti- 
cal lines were being measured and found that the position 
of the two was affected very little by the introduction of 
ellipticity. Therefore, we conclude that a spherical model 
is appropriate for our analysis. 

For every set of free parameters {M*/Ly, (3, 5 C , r sc }, 
we can compute the predicted position of the arcs, find the 
likelihood assuming gaussian distributions, and constrain 
the acceptable mass models. The largest radius at which 
the mass is probed is that corresponding to the location 
of the tangential arc (75.8 kpc). Now r sc is expected to 
be much greater than 100 kpc in CDM clusters (Bullock 
et al. 2001; see also Wu 2000). In this case, the location 
of the critical lines depends only marginally on r sc and 
the combined luminous and dark density profile has only 
3 free parameters (we fix r sc — 400 kpc in the following). 
Fig. 3 shows the likelihood contours (68%, 95% and 99%) 
in (3 vs. M*/Lv space, obtained using the likelihood ra- 
tio statistic after marginalization with respect to 8 C . Note 
that with lensing alone we can rule out with greater than 
95% confidence a M98 DM density profile {(3=1.5). 

3.3. Lensing + Dynamics 

The full power of our analysis is only realized when we 
combine the earlier constraints with those made by mea- 
suring the spatially-resolved stellar velocity dispersion pro- 
file of the BCG. Given our two-component mass model 
(Section 3.1), we solved the spherical Jeans equation (e.g. 
Binncy & Trcmaine 1987) assuming an isotropic velocity 
ellipsoid for the luminous component. 2 

The assumption of isotropy in the region probed (within 
30 % of R e ) we justify both on theoretical and observa- 
tional grounds. Numerical simulations (e. g. van Albada 
1982) and observations (e.g. Gerhard et al. 2001; Koop- 
mans & Treu 2002) appear to rule out significant tangen- 
tial anisotropy and permitting some radial anisotropy only 
at large radii. Strong radial anisotropy in the very central 
regions can also be ruled out based on the grounds of sta- 
bility (e. g. Mcrritt & Aguilar 1985; Stiavclli & Sparke 
1991) and consistency requirements (Ciotti 1999). How- 
ever, as a check, we ran Osipkov-Merritt (Osipkov 1979; 
Merritt 1985a, b) models with anisotropy radius set equal 
to progressively lower radii. Moving the anisotropy radius 
towards zero pushes the likelihood contours towards lower 
values of (3. 

For each set of parameters in the lens model, we com- 
puted the likelihood given the velocity dispersion profile, 
taking into account the effects of seeing, radial binning, 
and finite slit width. The total likelihood was computed 
by multiplying the velocity dispersion likelihood by the 

1 As a check of the results, we also considered a Hernquist (1990) luminous mass distribution. The results on /3 (see below) arc virtually 
unchanged, while slightly larger values of M/L for the stellar component are obtained. 

2 Spherical dynamical models have been shown to reproduce accurately the kinematics of slightly elongated galaxies like the BCG (e.g. 
Kronawitter et al. 2000). 



47rr 2 (r,/ + r) 2 ' 

mass density profile of total mass Mi, which reproduces 
well 1 the observed surface brightness profile (with R e = 
0.76r, 7 ). The DM halo is modeled as, 

Pc$c 



Pd(r) 



(2) 



(r/r sc y(l + (r/r sc jf-P> 

representing a generalization of the CDM-motivated ha- 
los, with an inner slope (3 (NFW and M98 correspond to 
(3 = 1,1.5 respectively; p c is the critical density). We 
assume that the BCG lies at the center of the overall po- 
tential. 

For a given stellar mass-to-light ratio M*/L v both M L 
and rj can be deduced from the surface photometry leav- 
ing 4 free parameters in our mass model: 1) M*/Ly] 2) 
the inner slope of the DM profile (3\ 3) the DM density 
scale S c ; and 4) the DM scale radius r sc . 

3.2. Gravitational Lensing 

Given our two-component spherical model, we adopted 
a simple lensing analysis using only the positions and red- 
shifts of the radial and tangential arcs (see Bartelmann 
1996). 

The locations of the radial and tangential arcs can be 
estimated by calculating the position of the corresponding 
radial and tangential critical curves of the projected mass 
distribution. The Jacobian matrix of the lens mapping has 
two eigenvalues, A r = 1 — 4- — and X t = 1 — I2 |, where 
x = R/r sc and m is a dimcnsionless function proportional 
to the mass inside projected dimensionless radius x (see, 
e. g., Bartelmann 1996, Schneider, Ehlers, Falco 1992). 
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likelihood obtained with the lensing analysis. The bot- 
tom panel in Fig. 3 shows the final results of the combined 
analysis. With 99% confidence, M^jLy lies between 2.3 
and 3.7 (broadly consistent with local values after passive 
evolution), with the inner DM slope (J5) lying between 0.05 
and 0.8, flatter than that expected from CDM simulations. 
The best-fitting parameters are M r /Ly = 3.1 M Q /L y, 
(3 = 0.35, 5 C = 24000. 
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Fig. 3.— Likelihood contours (68%, 95%, and 99%) obtained 
for the mass modeling of MS2137-23 with a Jaffe luminous distri- 
bution plus a generalized NFW DM distribution. (Top): Contours 
obtained from the position of the radial and tangential arcs alone. 
Note that a M98 (/3 = 1.5) profile is excluded at the 95% level. (Bot- 
tom): Contours obtained including the measured velocity dispersion 
profile. Note the improved constraints on the mass parameters and 
that NFW profiles are clearly ruled out at the 99% level. 

To check our results, we changed R e by 10% with a 
negligible effect on the likelihood contours. Similarly, the 
contours are virtually unchanged by offsets of C/5 in the 
position of the tangential arc and by changing the seeing 
by 30%. Changing the position of the radial critical line 
by ±0'/5 shifts the likelihood contours by ±0.1 in the (3 
direction. 

Systematic offsets of the velocity dispersion profile due 
to template mismatch and poor continuum fitting intro- 
duce correlation between the kinematic points that are not 
considered in the likelihood ratio analysis. To investigate 
this, we have repeated the analysis with the velocity dis- 
persion profile shifted by the estimated systematic error 
(~ 15 km s _1 ). A lower overall velocity dispersion profiles 
shifts our likelihood contours towards lower (3 {[3 < 0.65 
99 % CL), and viceversa ((3 < 0.9 99 %CL): our hard 99% 



CL upper limit is f3 < 0.9. 

4. SUMMARY AND DISCUSSION 

We have secured spectroscopic redshifts for the radial 
and tangential arcs in the cluster MS 2137-23 and deter- 
mined the velocity dispersion profile of the BCG to a phys- 
ical radius of ~8 kpc. We have combined these measure- 
ments with the lensing geometry in order to construct a 
self-consistent model of the mass distribution in the clus- 
ter core. Using a spherically symmetric luminous and DM 
mass distribution, we rule out the presence of a DM halo 
with an inner slope (3 > 0.9 at greater than 99% confi- 
dence, including systematics. M98 and NFW-type halos 
with (3 > 1.0 are inconsistent with the mass distribution 
in the core of MS2137-23. 

Since the infall of baryons associated with the BCG are 
likely to steepen the DM halo (Blumenthal et al 1986; Mo, 
Mao and White 1998), our measured profile may imply 
the original DM profile was even flatter. A full model- 
ing of this process is beyond the scope of this letter, but 
this strengthens our conclusion that the inner regions of 
the DM halo of MS2137-23 cannot be described by CDM- 
motivated universal halos. A potential concern is that our 
models have only two mass components, but X-ray emit- 
ting gas could be a non-negligible third massive compo- 
nent. Using ROSAT observations of MS2137-23 (Ettori & 
Fabian 1999) we estimate that removing the X-ray com- 
ponent will steepen the resulting DM halo slope by less 
than ~ 0.1 and therefore does not change dramatically 
the result. 

Finally, individual halo shapes can depart from the en- 
semble average behavior. Therefore it is necessary to apply 
such a test to a sample of clusters. Our simple method is 
applicable to all approximately round clusters with a mas- 
sive galaxy at their center, provided that they have at 
least a giant tangential arc (radial arcs further enhance 
the sensitivity but are not required). We are in the pro- 
cess of collecting data for a dozen clusters with the aim of 
performing such a statistical test. 
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